Kinematic Evidence of Minor Mergers in Normal Sa Galaxies: 
NGC 3626, NGC 3900, NGC 4772 and NGC 5854 
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ABSTRACT 

BVRI and Ha imaging and long-sht optical spectroscopic data are presented for 
four morphologically normal and relatively isolated Sa galaxies, NGC 3626, NGC 3900, 
NGC 4772 and NGC 5854. VLA H I synthesis imaging is presented for the first three 
objects. In all four galaxies, evidence of kinematic decoupling of ionized gas compo- 
nents is found in the long-slit spectroscopic data; the degree and circumstances of the 
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distinct kinematics vary from complete counterrotation of all of the gas from all of the 
stars (NGC 3626) to nuclear gas disks decoupled from the stars (NGC 5854) to anoma- 
lous velocity central gas components (NGC 3900 and NGC 4772). In the three objects 
mapped in H I, the neutral gas extends far beyond the optical radius, RH1/R25 > 2. In 
general, the H I surface density is very low and the outer H I is patchy and asymmetric 
(NGC 3900) or found in a distinct ring, exterior to the optical edge (NGC 3626 and 
NGC 4772). While the overall H I velocity fields are dominated by circular motions, 
strong warps are suggested in the outer regions by bending of the minor axis isoveloc- 
ity contours (NGC 3900) and/or systematic shifts in position angle between inner and 
outer rings (NGC 3626 and NGC 4772). In the interior, coincidence is found between the 
Ha and H I in rings, sometimes partial and crisscrossed by dustlanes. Optical imaging 
is also presented for NGC 4138 previously reported by Jore et al. (1996) to show coun- 
terrotating stellar components. The multiwavelength evidence is interpreted in terms 
of the kinematic "memory" of past minor mergers in objects that otherwise exhibit no 
morphological signs of interaction. 

Key words: - galaxies: evolution - galaxies: individual (NGC 3626, NGC 3900, 
NGC 4138, NGC 4772, NGC 5854) - galaxies: interactions - galaxies: kinematics and 
dynamics - galaxies: spiral 
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1. Introduction 

While mergers between comparable mass galaxies 
may be responsible for some of the most dramatic 
extragalactic events, minor mergers with small mass 
satellites may play a subtle but nonetheless critical 
role in the evolution of garden- variety spirals, in- 
cluding perhaps the Milky Way. It is easy to rec- 
ognize the fireworks of the former by the tides, star- 
bursts and nuclear activity they trigger. The latter 
however may lurk hidden to morphological inspec- 
tion. Minor mergers are suggested to produce the 
X-shaped structures seen in some peculiar SO's (Mi- 
hos et al. 1995), the departures from axisymmetry 
seen in a significant fraction of spiral disks (Zaritsky 
& Rix 1997), the driving mechanism for strong stellar 
bars (Laine & Heller 1999) and the peculiar extended 
counter-rotating disks seen in some early -type spi- 
rals, e.g. NGC 7217 (Merrifield & Kuijken 1994; Buta 
et al. 1995), NGC 3593 (Bertola et al. 1996; Corsini et 
al. 1998), NGC 3626 (Ciri et al. 1995) and NGC 4138 
(Jore et al. 1996). 

Satellites are commonly found in the vicinity of 
normal galaxies (Zaritsky et al. 1997). If the Local 
Group is at all typical of spiral extragalactic environ- 
ments, then the characteristics of such satellites are 
familiar. Interaction of the Milky Way with its low 
mass satellites is hypothesized to explain a variety of 
local kinematic peculiarities: the thick disk, the outer 
flare and warp, the Magellanic Stream and the high 
velocity H I clouds. 

Several authors have begun to address numerically 
the complex phenomenon of minor mergers of gas- 
rich satellites with more massive spiral companions. 
An immediate worry is how to avoid over-heating the 
inner disk of the primary or disrupting it entirely. 
Minor mergers of low-density, small mass, gas rich 
companions may avoid vertical heating, thereby pre- 
serving the victim disk. Predictions of the observable 
consequences of the accretion of a small mass compan- 
ion by a spiral include the warping or actual spreading 
out of the disk (Quinn et al. 1993). As the satellite 
material sinks towards the center of the large galaxy, 
the gas, unlike the stars, may lose angular momentum 
because of dissipation (Hernquist & Mihos 1995). The 
resultant inflow of gas towards the center may trig- 
ger a nuclear starburst (Mihos & Hernquist 1994) or 
fuel nuclear activity (Taniguchi 1999). Because of the 
predicted heating, even a satellite with one tenth the 
mass will thicken the disk and possibly build-up the 



bulge, driving the morphology of the post merger pri- 
mary towards an earlier spiral type. The discovery 
of kinematically decoupled disks in early spirals may 
thus prove the efflcacy of minor merger models. 

Intrinsic galaxy properties arise from two principal 
components: absolute "scale" (size, luminosity, mass) 
and "form" (morphology, bulge-to-disk ratio, color, 
gas content, star formation rate). It is well-known 
(e.g. Roberts & Haynes 1994) that the early spiral 
types show a much greater spread in their character- 
istic properties related to "form" than do their later- 
type counterparts. As a part of a study to investigate 
the heterogeneity and dark matter content of the Sa 
galaxy class, we have investigated the kinematics of 
gas and stars in a sample of morphologically normal, 
isolated galaxies classified as Sa in the Revised Shap- 
ley Ames Catalog (Sandage & Tammann 1987; RSA). 
In the course of that study, we have obtained major 
and minor axis long slit spectra for both stars and 
gas for a sample of 20 objects; for nine of them, we 
have also obtained H I synthesis maps yielding also 
two-dimensional H I velocity fields. A summary of 
results focussing on the data presentation and mass 
modelling was presented in Jore (1997). The analy- 
sis of the mass modelling and stellar velocities will be 
presented elsewhere (Jore et al., in preparation). 

During the course of that work, a number of un- 
usual cases of kinematic decoupling were discovered. 
In Jore et al. (1996), we discussed the distinctive 
case of NGC 4138, an isolated Sa with two extended 
coplanar counter-rotating stellar disks embedded in 
a huge HI disk that co-rotates with the secondary 
stellar component. In a separate paper (Jore et 
al. 2000), we discuss five Sa galaxies, NGC 1169, 
3623, 4866, 5377 and 5448, which represent a range of 
environments and morphologies within the Sa class. 
Of those objects, NGC 1169 and NGC 5448, ex- 
hibit no pecTiliar optical kinematics but their outer 
H I disks are strongly warped and kinematically 
skewed. NGC 3623, whose companions in the Leo 
Triplet are known to be interacting, exhibits a trun- 
cated H I disk and central gas decoupling. NGC 4866, 
viewed nearly edge-on, and NGC 5377 show dou- 
bling of their central gas components, as well as outer 
H I warps. Here, we present four additional examples: 
NGC 3626, NGC 3900, NGC 4772 and NGC 5854. 
Along with NGC 4138, each of these morphologically 
normal Sa galaxies exhibits peculiar kinematics that 
betray a disturbed past. We propose that each rep- 
resents a different stage and circumstance of involve- 
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ment in a relatively minor merger event with a gas- 
rich companion that is no longer recognizable. 

This paper presents evidence of this past interac- 
tion through a combination of optical imaging and 
spectroscopy and H I synthesis mapping. In Section 2, 
we discuss the acquisition and reduction of the body 
of data: broad- and narrow-band optical images, op- 
tical long-slit spectra, and H I synthesis maps. The 
kinematic evidence of past interaction is presented for 
each galaxy in Section 3. Section 4 gives a discussion 
of the observational results in the context of the minor 
merger scenario. 

2. Observations 

The original observing program was designed to 
provide the details of the galaxy dynamics for the 
purpose of modelling the dark matter distribution in 
Sa galaxies and, as such, consisted of long-slit major 
and minor axis spectra of both stars and ionized gas. 
For a subset, principally gas-rich objects showing evi- 
dence in the optical spectra of kinematic peculiarities, 
H I synthesis mapping and photometric broad- and 
narrow- band imaging was also undertaken. 

The basic properties of these four galaxies plus 
NGC 4138 discussed by Jore et al. (1996) are summa- 
rized in Table 1. The first entries include the (B1950) 
R.A. and Dec. and the morphological classifications 
as given in both the RSA (Sandage & Tammann 1987) 
and the Third Reference Catalog of Bright Galaxies 
(de Vaucouleurs et al. 1991; RC3). A morphological 
description of the Ha emission is then given, followed 
by a set of radial measures: D25 x ^25, in arcminutes, 
from the RC3 and the corresponding radius R25, in 
arcseconds; Rb25, in arcseconds, derived from the B 
image; R/83L, in arcseconds, the radius encompass- 
ing 83% of the total / band light, following Haynes 
et al. (1999a); and Rorc, in arcseconds, the maxi- 
mum radius at which the optical rotation curve is 
derived, for any species, stars or gas. The corrected 
integrated blue magnitude, BSJ,, and the color indices, 
(U—B) and {B — V), all from the RC3, are then listed, 
along with the optical heliocentric velocity, Vq, in 
km s~^, and the adopted distance, in Mpc. Since all of 
these galaxies are nearby, distances are adopted from 
available estimates accounting for group membership 
and local peculiar motions and/or supernova distance 
moduli. The total B luminosity, Lb, is derived from 
Bj, and the distance. The total apparent magnitudes 
ms, my, m^^, and ni/, and their estimated errors. 



have been derived from the images following Haynes 
et al. (1999a). Likewise, the ellipticity e = (1 — ^) and 
position angle are the mean disk properties extracted 
from the / band isophotal fits using the "marking the 
disk" method of Haynes et al. (1999a). Finally, the 
adopted R band bulge-to-disk ratio B /D/j is given. 

The acquisition and reduction of the multiwave- 
length data sets are discussed here separately. 

2.1. Optical Images 

BVR and Ha images were obtained with the KPNO 
0.9m telescope in its f/7.5 configuration over the 
nights of April 19-23, 1996. Additionally, / band im- 
ages, some taken under non-photometric conditions, 
were obtained with the same instrument by R. Gio- 
vanelli and MPH as part of another project. For all 
images, the 2048x2048 pixel T2KA CCD was used. 
Its read noise is 4 e~ and the gain was set to 5.4 e~ 
ADU~^, except for the the / band images, when the 
gain was set to 3.6 e~ ADU~^. With the field correc- 
tor in, the images had a pixel scale of 0.68". In this 
section, we describe the image acquisition, processing 
and analysis techniques. 

2.1.1. Broadband BVRI Images 

The broadband images were exposed over the whole 
chip, giving a total field of view of 23' on a side. 
The images were trimmed, bias subtracted, and flat 
fielded using CCDPROC in IRAF 0. Flat fields for 
_B, V and R were obtained from twilight sky observa- 
tions; for the / band, a fiat field was constructed from 
the median of 30-50 object frames taken throughout 
the night. Further image reductions made use of the 
GALPH0T| package of IRAF/STSDAS scripts for 
galaxy surface photometry (see Haynes et al. 1999a). 
Regions of the image free of galaxy emission and 
bright stars were interactively marked. An average 
value of the sky contribution in these boxes after 
masking of the stars within them was subtracted as 
the sky contribution. Images were also edited inter- 
actively to remove stars and cosmic rays and con- 
volved to guarantee the same resolution before com- 
bining. Images with three or more observations on 

^IRAF is distributed by the National Optical Astronomy 
Observatories. 

^The GALPHOT surface photometry package is a collection 
of IRAF/STSDAS scripts originally developed by Wolfram 
Freudling and John Salzer; the Cornell version has been further 
modified and is maintained by MPH. 
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Table 1 
Galaxy Properties 





IVr QROR 
IN ODZD 


IN oyuu 


M /II QQ 
IN ^100 


M /I77O 
IN i Zi 


IVr KQK/I 
IN OoO^ 




TT ^^Q/t Q 

U Do4o 


TT ^TQ^ 

U blab 




TT cnoi 


TT 070^^ 
U y / ZD 


R.A.(B1950) 


111725.9 


114633.3 


120659.3 


125055.9 


150516.2 


Dec.(B1950) 


+183756 


+271806 


+435757 


+022627 


+024537 


RSA Type'' 


Sa 


S(r)a 


S(r)a pec 


Sa: 


Sa 


RC3 Type'' 


RLAT+.. 


.LAR+.. 


.LAR+.. 


.SASl.. 


LBS+./ 


Ho; morph 


n+cr 


(trace) 


n+cb+or 


n+cb+or 






2.69x 1.95 


3.16 X 1.70 


2.57 X 1.70 


3.39 X 1.70 


2.75 X 0.79 




81 


96 


78 


102 


84 


I^-B25 ( ) 


76 


97 


85 


132 


78 


R/83L (") 


56 


75 


54 


120 


48 


^^orc \ } 


45 


53 


50 


58 


51 




11.75 


12.24 


12.14 


11.90 


12.68 


{U - Bf 


0.30 




0.31 


0.26 


0.29 


(B - V)^ 


0.81 


0.82 


0.83 


0.83 


0.81 


V© (km s"^)"" 


1493 


1798 


888 


1040 


1663 


Distance (Mpc) 


24 


27 


16 


16 


24 


Lb (IQio Lq) 


1.8 


1.4 


0.6 


0.7 


0.8 


ms 


11.91 ±.04 


12.40 ±.02 


12.12 ±.04 


11.97 ±.04 




my 


11.11 ±.02 


11.65 ±.02 


11.31 ±.02 


11.15 ±.03 


11.20 ±.10 




10.59 ±.02 


11.08 ±.02 


10.73 ±.03 


10.58 ±.02 


10.64 ±.02 


m/ 




10.42±.02 


10.03±.02 






e 


0.32±0.02 


0.48±0.02 


0.39±0.02 


0.48±0.03 


0.67±0.02 


RA. (°) 


164±3 


4±1 


150±1 


148±2 


55±1 


B/Dfl 


0.37 


0.27 


0.19 


0.44 


0.82 



''From RSA. 
''From RC3. 

'^Code for morphology: n=miclcar; cr=corc ring; cb=elongated core; or=outer ring. 
''Maximum extent of optical rotation curve traced by any species. 
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the same filter were scaled and median combined to 
remove cosmic rays. For images with fewer exposures 
cosmic rays were removed before combining using the 
task COSMIC in IRAF and interactively editing all 
remaining cosmic rays by hand. The filters, total ex- 
posure times, and the seeing values for the combined 
images are presented in Table 2 for each galaxy. 

Elliptical isophotal fitting routines were run on the 
star edited BVRI images of each galaxy. These rou- 
tines require the user to mark an initial major and 
minor axis. The program then fits isophotal ellipses 
moving outwards and inwards from the marked initial 
location and integrates the flux within the ellipses. 
Using these results, disk parameters such as position 
angle, ellipticity, and slope are fit and extrapolated 
magnitudes are found. Selected values of these pa- 
rameters are included in Table 1. 

Conversion from instrumental to apparent magni- 
tudes was accomplished using frequent observations 
of standard stars (Landolt 1992) observed throughout 
each night. The photometric coefficients for two pho- 
tometric nights during which the BVR images were 
obtained (nights 1 and 3) are listed in Table 3. Pho- 
tometric coefficients were also determined for night 
2; however, they fail to reproduce Landolt standard 
magnitudes accurately and were therefore discarded 
as non-photometric. All images from night 2 were 
subsequently scaled to ones obtained on nights 1 or 3 
in order to insure photometric accuracy. 

Color plots of each galaxy in B — V, B — R and 
V — R were produced by first aligning all the band- 
passes and convolving them to the same resolution 
(the poorest of the set, as given in Table 2). The 
images were then converted into instrumental mag- 
nitudes and rcscalcd to apparent magnitudes using 
the appropriate photometric coefficients. Final color 
images were made by subtracting the resultant im- 
ages and correcting for galactic extinction. Values 
for galactic extinction were calculated using A(,(B) 
values from the RC3 and the conversions to other 
bands given in Ricke & Lebofsky (1985). Addition- 
ally, model color images were also created from the 
broadband isophotal ellipse fitting routines mentioned 
above. Using the results of the ellipse fitting, model 
broadband images were created, converted into ap- 
parent magnitudes as above, and subtracted to form 
noiseless model B — V, B — R and V — R images. 

The separation of bulge and disk components was 
accomplished using FITMAP, a two-dimensional de- 
composition routine written and kindly provided by 



G. Moriondo (Moriondo et al. 1998). FITMAP fits a 
model with a specified bulge and exponential disk to 
the data, adjusting the allowed parameters until the 
optimum model is found as judged by the reduced 
value. The bulge shape index n for the generalized 
exponential law Ib{R) = /eea;p[— (i?/i?e)^^"] can be 
specified as an integer value, with n = 1, an exponen- 
tial bulge and n = 4, a de Vaucouleurs bulge. In order 
to speed processing time, the fit was performed only 
on an averaged quadrant of the galaxy. In this data 
set, the average was typically a combination of all four 
quadrants, achieved by folding the image along both 
major and minor axes, except where bright stars were 
found close to the galaxy center. All filters were fitted 
separately to allow comparison of the results obtained 
with different bandpasses. 

Recent literature (e.g. de Jong 1996, Moriondo 
et al. 1998) shows that for earlier type spirals, such 
as those used in this sample, an n = 2 bulge is the 
best fit. That finding is supported likewise here. Fur- 
thermore, the results indicate an increase in the B /D 
parameter with redder wavelengths, compatible with 
the bulge having an overall redder color than the disk 
and there being less obscuration from dust at longer 
wavelengths. A summary of the B/D parameters ob- 
tained for each filter and adopting the n = 2 bulge is 
given in Table 4. 

2.1.2. Narrowband Ha Im,ages 

Ha images were produced by observing the same 
region in two narrow (SOA) bandpasses, one centered 
on the appropriately redshifted Ha emission line (the 
ON image) and a second just off the emission line 
(the OFF image). Both OFF and ON images were 
obtained over the nights of April 22-23, 1996 with 
the KPNO 0.9m telescope. The same T2KA chip, as 
described above, was used; however, only the central 
1024x1024 pixel region was exposed. Usually, two ON 
exposures of 30 minutes each and a single OFF one 
of the same duration were obtained; Table 2 gives 
the details of the Ha observations for each galaxy. 
Images were trimmed, bias subtracted, fiat fielded, 
and cosmic ray edited in the same manner as the BVR 
images described above. The ON and OFF images 
were then convolved to the same resolution, and the 
two ON images were combined. The OFF image was 
scaled to the summed ON image and subtracted to 
produce an Ha line image. The nights were both 
non-photometric and no calibration is available. 



6 



Table 2 

Details of Imaging Observations 



Object 


Filter 


Night 


(sec) 


Seeing 
tt 


Note 


N3626 


B 


1 


3600 


2.69 


a 




V 


1 


1350 


2.40 


a 




R 


1 


1080 


1.97 






I 




600 


2.99 


b 




ON 


4 


3600 


1.98 


b 




OFF 


4 


1800 


1.98 


b 


N3900 


B 


3 


3600 


2.03 






V 


3 


1350 


2.02 






R 


3 


900 


2.16 






I 




600 


1.62 






ON 


5 


3600 


1.81 


b 




OFF 


4 


1800 


1.77 


b 


N4138 


B 


1 


3600 


2.24 






V 


1 


1350 


1.86 






R 


1 


900 


1.81 






I 




600 


2.12 






ON 


5 


3600 


2.00 


b 




OFF 


4 


1800 


2.09 


b 


N4772 


B 


1 


3600 


1.57 






V 


3 


1350 


1.65 






R 


3 


900 


1.65 






I 




600 


1.92 


b 




ON 


5 


3600 


1.56 


b 




OFF 


5 


1800 


1.64 


b 


N5854 


B 


1 


3600 


2.50 


b 




V 


1 


1350 


2.19 






R 


1 


900 


1.99 






I 




600 


2.58 


b 



''Tracking problems; image quality poor 
Nonphotometric 
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Table 3 

Photometric Solutions for BVR Imaging 



Night 


Filter 


Zoropoint 


Color Coef. 


Extinc. Coef. 






(mag) 






1 


B 


21.005 (0.039) 


0.099 (0.019) 


0.237 (0.004) 




V 


21.311 (0.018) 


0.000 (0.009) 


0.157 (0.004) 




R 


21.411 (0.013) 


0.009 (0.014) 


0.104 (0.003) 


3 


B 


21.062 (0.007) 


0.075 (0.012) 


0.233 (0.004) 




V 


21.354 (0.013) 


-0.017 (0.010) 


0.155 (0.003) 




R 


21.424 (0.010) 


0.019 (0.016) 


0.102 (0.003) 



Table 4 

Results of 2D Bulge-Disk Decomposition. 







Exponential Disk 


n = 2 Bulge 








Object 


Filter 


(mag / arcsec?) 


h 

(arcsec) 


{mag / arcsec?) 


(arcsec) 




B/D 




N3626 


B 


20.94(0.01) 


21.3(0.16) 


18.94(0.05) 


3.2(0.08) 


2.73 


0.37 






V 


20.10(0.01) 


21.1(0.12) 


17.99(0.05) 


3.0(0.06) 


3.75 


0.37 






R 


19.55(0.01) 


20.7(0.13) 


17.50(0.04) 


3.0(0.06) 


3.20 


0.37 






I 




21.3(0.34) 




3.4(0.05) 


1.16 


0.30 


a 


N3900 


B 


21.68(0.02) 


24.4(0.21) 


20.69(0.05) 


4.8(0.18) 


1.55 


0.25 






V 


20.78(0.01) 


23.0(0.16) 


19.77(0.04) 


4.9(0.13) 


1.61 


0.30 






R 


20.15(0.01) 


22.2(0.14) 


19.16(0.04) 


4.8(0.12) 


1.84 


0.27 






I 


19.27(0.01) 


20.7(0.11) 


18.25(0.02) 


4.0(0.06) 


0.79 


0.25 




N4138 


B 


20.46(0.03) 


16.6(0.14) 


21.44(0.08) 


7.0(0.48) 


2.19 


0.19 






V 


19.67(0.01) 


16.9(0.08) 


19.90(0.04) 


5.2(0.16) 


1.64 


0.20 






R 


19.06(0.01) 


16.7(0.06) 


19.01(0.04) 


4.4(0.11) 


1.72 


0.19 






I 


18.37(0.01) 


17.0(0.05) 


17.95(0.02) 


3.7(0.06) 


1.02 


0.18 




N4772 


B 


22.81(0.01) 


50.4(0.36) 


21.22(0.01) 


9.6(0.10) 


0.44 


0.41 






V 


21.96(0.01) 


48.4(0.25) 


20.27(0.01) 


9.3(0.08) 


0.65 


0.45 






R 


21.35(0.01) 


46.8(0.23) 


19.64(0.01) 


9.1(0.07) 


0.72 


0.47 






I 




42.3(0.40) 




8.1(0.07) 


0.23 


0.44 


a 


N5854 


B 




20.0(0.38) 




6.0(0.14) 


0.66 


0.63 


a 




V 


20.15(0.22) 


19.5(0.35) 


19.07(0.23) 


5.7(0.12) 


0.55 


0.61 


a 




R 


19.59(0.02) 


19.5(0.15) 


18.54(0.03) 


5.7(0.12) 


0.73 


0.59 






I 




19.6(0.17) 




5.6(0.11) 


1.05 


0.63 


a 



''Nonphotometric 
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2.2. Optical Long Slit Spectra 

These data are described in more detail in Jore 
(1997). Optical long slit spectra were taken at the 
Palomar 5m Hale Telescope^ using the red and blue 
cameras of the Double Spectrograph. The slit dimen- 
sions were 120" by 2". A dichroic at 5500 A and grat- 
ings in both cameras with dispersions of 1200 lines 
mm~^ were used. The red camera was centered on 
the emission lines of Ha and [N II] A 6548 A and A 
6584 A, with the [S II] A 6717 A and A 6731 A lines 
also present. The blue camera was centered on the 
Mg I b A 5183.6, 5172.7 and 5167.3 A stellar absorp- 
tion lines with the [O III] A 5007 A emission line also 
visible. Table 5 summarizes the Double Spectrograph 
setup, while Table 6 presents the details of all spectral 
observations used to extract rotation curves. All ex- 
posures over 40 minutes composed of co-added shorter 
exposures. It should be noted that some spectra were 
taken under conditions of poor transparency so that 
the effective integration is significantly less than the 
actual exposure time. Arc lamps and radial velocity 
standard stars were observed for wavelength and ve- 
locity calibration. In general, arc lamps were used for 
wavelength calibration on the blue side, while night 
sky lines provided the necessary scaling anchor on the 
red side. 

Rotation curves were extracted for each of the ion- 
ized gas components using a Gaussian fitting rou- 
tine. The stellar rotation curves were extracted us- 
ing several different methods employing both cross- 
correlation and direct-fitting techniques. Further dis- 
cussion of the reduction procedure is presented in Jore 
et al. (1996) and Jore (1997). The gas rotation curves 
displayed in Section 3 are derived from Ha and/or 
[Nil] , depending on the relative strengths of the lines 
and the possible confusion with line absorption in the 
case of Ha. The stellar rotation curves are discussed 
here; the stellar velocity dispersions will be discussed 
elsewhere (Jore et al., in preparation). 

2.3. HI Synthesis Imaging 

The H I 21~cm line observations were conducted 
with the Very Large Array (VLA)^ in its C configu- 
ration. As discussed in Section 3.4, NGC 5854 is very 



Observations at the Palomar Observatory were made as part 
of a continuing collaborative agreement between the California 
Institute of Technology and Cornell University. 
'The Very Large Array is a facility of the National Radio As- 
tronomy Observatory. 



gas-poor and was not mapped. A summary of the 
observational setup for the VLA observations of the 
other three galaxies is given in Table 7. Each galaxy 
was allocated 8 hours of array time in its C configura- 
tion, with a nominal beam size of ^--^12" and sensitivity 
to structures smaller than 7'. Within each 8 hour al- 
location, a flux calibrator was observed three times 
for ^10 minutes, and a phase calibrator was observed 
for 4 minutes between each 50 minute observation on 
the galaxy. All observations were obtained using the 
4ABCD spectral line mode with two overlapping spec- 
tral bands of 32 channels per band, with an overlap 
of at least 10 channels so that the total number of 
channels was 54, each 20.6 km s^^ wide. The total 
amount of on-source integration time and the total 
velocity coverage of the combined spectral bands are 
listed in Table 7. 

The UV data were flux and phase calibrated us- 
ing the AIPS software package after flagging bad vis- 
ibilities. The bandpass for each antenna of the line 
data was calibrated for each IF separately. The two 
spectral bands were then combined using the task 
UVGLU. Dirty maps were made using MX and were 
inspected to determine the channels free of line emis- 
sion. The continuum was subtracted in the UV plane 
using UVLIN by interpolating the continuum from 
the emission-line free channels to those with line emis- 
sion. Cleaned map cubes were created from the UV 
data using a Gaussian beam derived from the FWHM 
of the dirty beam. For most of the galaxies, two maps 
with different resolutions were made; one (using natu- 
ral weighting) with lower spatial resolution and higher 
signal-to-noise, and the other (using uniform weight- 
ing) with higher resolution but noisier. The beam 
size and position angle and RMS noise per channel 
for each map are listed in Table 7. 

Once the map cubes were made, further analysis 
of the data was done using the GIPSY (van der Hulst 
et al. 1992) data reduction package. The data cubes 
were first smoothed to twice the original beam size 
and then clipped at the la level. The regions contain- 
ing line emission in two or more consecutive channels 
were interactively extracted. The data in the original 
cube corresponding to the regions extracted in the 
smoothed cube were conditionally transferred. This 
final cube contains only the line emission from the 
sources. Zeroth and first order moment maps (inte- 
grated intensity and velocity field maps) were made 
using the GIPSY task MOMENTS. 

Column density maps were obtained from the in- 
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Table 5 

Double Spectrograph Observing Setup 



Dates 


Camera 


Chip 


Spat. Res. 


A Range 


A Res. 


Vel. Res. 








(" Pix-i) 


(A) 


(A pix-i) 


(km s-i) 


Mar. 94 - 


Red 


TI 167 


0.58 


6190-6850 


0.82 


38 


Apr. 95 


Blue 


TI 432 


0.78 


4950-5400 


0.56 


38 


Mar. 96 


Red 


1396CR14-0 


0.46 


6150-6825 


0.66 


30 




Blue 


TI 432 


0.78 


4950-5400 


0.56 


38 



Table 6 

Optical Long-Slit Spectral Observations 



Object 


Date 


Total Int. 
(seconds) 


PA 

(° E of N) 


Axis 


NGC 3626 


April 1995 


3600 


157 


major 






1800 


67 


minor 


NGC 3900 


March 1994 


3600 


2 


major 






3600 


92 


minor 




April 1995 


3600 


2 


major 






1800 


92 


minor 


NGC 4772 


April 1995 


3600 


147 


major 






1800 


57 


minor 


NGC 5854 


March 1994 


3600 


55 


major 






3000 


145 


minor 




April 1995 


3600 


55 


major 






1800 


145 


minor 




March 1996 


6419 


55 


major 
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Table 7 

VLA C HI Synthesis Observing Parameters 





weighting 


N 3626 


N 3900 


N 4772 




scheme 


U 6343 


U 6786 


U 8021 


Date Obs. 




2/23/96 


2/24/96 


11/10/94 


Tint (min) 




373.5 


371.5 


403 


Velocity Coverage (km s""') 




930 - 2012 


1251 - 2336 


485 - 1564 


Beam Size (") 


uniform 


12.1 X 11.8 


11.7 X 11.7 


13.8 X 12.4 




natural 


16.8 X 16.2 


22.5 X 20.0 


31.7 X 23.1 


Beam PA (°) 


uniform 


113 


165 


143 




natural 


101 


79 


137 


RMS channel"^ (mjy beam~^) 


uniform 


0.28 


0.36 


0.44 




natural 


0.23 


0.31 


0.33 


minimum Nhi (10^" atoms cm~^) 


natural 


1.0 


0.3 


0.2 



tegrated intensity maps using the relation 

N{H) = 1.10297xlO^\l+v/cf^j S{v)dv [atoms cm-^] 

(1) 

where J S{v)dv is the integrated flux (in Jy/beam) 
times the velocity width (in km /channel), a and 
b are the FWHM of the major and minor axes of the 
beam (in arcscc) , and v is the systemic velocity of the 
source (in km s~^). The minimum detectable column 
density quoted in Table 7 is the 2a value in the natural 
weight column density map. 

2.4. Continuum emission 

Cleaned maps were made of the non-continuum 
subtracted maps. The continuum map for each source 
was determined by taking the mean of all line-free 
channels. The primary beam corrected flux density 
of nearby sources was calculated using the GIPSY 
routine FLUX. The continuum results for each galaxy 
are discussed separately: the derived continuum flux 
density is included in Table 8. 

2.4.1. NGC 3626 

For NGC 3626, the average of 17 line-free channels 
was calculated, yielding a continuum map. Contin- 
uum emission was detected at the center of NGC 3626 
with a flux density 5*1413 — 9.9 mJy, slightly larger 
than the flux density 6*1415 = 6 mJy measured by 
van Driel et al. (1989). There is another nearby 
source located 2^9 west of the center (ll''17™13'!5, 



+18°37'38" (1950)) with a flux density 6*1413 = 4.1 mJy. 
The companion galaxy, UGC 6341, has no continuum 
emission associated with it, at a 3(t point source sen- 
sitivity of 0.3 mJy. 

2.4.2. NGC 3900 

For NGC 3900, the continuum map was constructed 

similarly. Continuum emission was detected at the 
center of NGC 3900 with a flux density 61412 = 
1.2 mJy in agreement with the upper limit of 6*1415 < 
6 mJy given by van Driel et al. (1989). A nearby 
source is located 2.'8 north of the center (ll''46'"35!3, 
+27°20'24" (1950)) with a flux density 61412 = 1-4 mJy 
which is coincident with a small background galaxy 
seen in the optical images. 

2.4.3. NGC 4772 

For NGC 4772, the average of 16 line-free channels 
was used to produce the continuum map. Continuum 
emission was detected at the center of NGC 4772 with 
a flux density 61415 = 3.4 mJy. Four other nearby 
sources were found. Their positions and flux densities 
are as follows: (a): 12^50"59!9, +02°24'43" (1950) 
with a flux density 61415 = 5.0 mJy; (b): 12''51™06!0, 
-F02°26'39" (1950) with a flux density 61415 = 4.0 mJy; 
(c): 12''51"09!5, +02°28'15" (1950) with a flux den- 
sity 61415 = 2.3 mJy; and (d): 12''50'"58!3, +02°29'35" (1950) 
with a flux density 61415 = 3.7 mJy. 
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2.5. HI Channel Maps and Global Profiles 

The channels maps extracted from the natural 
weight data cubes are shown in Figures 1 to 3 for 
the individual galaxies NGC 3626, NGC 3900 and 
NGC 4772 respectively. The global HI profiles ex- 
tracted for the three target galaxies and for two oth- 
ers, UGC 6341 and NPMlG-l-18.0219, both found in 
the field of NGC 3626, are shown in Figure 4. Table 8 
presents the properties of these five objects as derived 
from the VLA maps; NGC 4138 is also included for 
comparative purposes. Entries are: (1) the contimmm 
flux density derived from the VLA map as discussed in 
Section 2.4, in mJy; (2) the H I deficiency, < DEF >, 
given by Magri (1994). < DEF > is the logarithmic 
difference between the observed H I mass and that 
expected for an isolated galaxy of the same morpho- 
logical class and linear optical diameter (Haynes & 
Giovanelli 1984); (3) the H I line flux / Shi dV in 
Jy-km s~^. For the VLA datasets, / Shi dV is de- 
rived from the natural weight data cube, after appli- 
cation of a primary beam correction. (4) the H I mass, 
derived from / Shi dV, in M©; (5) the H I mass to 
blue luminosity ratio, M hi I Lb, in solar imits; (6) the 
H I line heliocentric systemic velocity, V21, in km s~^; 
(7) the width of the H I line global profile at 20% 
of the peak intensity, AV21; in km s~^; (8) the ratio 
of the H I to optical radius, RH1/R25, where Rhi is 
measured at a constant surface density of 1 Mq pc~^; 
(9) the position of the dynamical center; (10) the incli- 
nation derived from the tilted ring model, in degrees; 
(11) the maximum rotational velocity derived from 
the best-fit rotation curve, Vmax, in km s~^. Further 
discussion is deferred to the next section in which each 
galaxy is discussed individually. 

3. Results for Individual Objects 
3.1. NGC 3626 

NGC 3626 is found in an extended cloud of galax- 
ies related to the NGC 3607 group, LGG 237 (Garcia 
1993); the central concentration around NGC 3607 is 
located about 50' (350 kpc at the adopted distance of 
24 Mpc) away. The only object known to be associ- 
ated with the group that lies within 30' of NGC 3626 
is UGC 6341, a Magellanic spiral at similar redshift, 
located at a separation of 5^9, 41 kpc. UGC 6341 
was also detected in the VLA H I synthesis map and 
is discussed further below. A ten minute exposure 
with the Palomar 5m telescope and the Double Spec- 
trograph of the faint galaxy located at R. A. (1950) = 



ll''17"'30!l, Dcc.(1950) = +18°34'31" failed to reveal 
any lines and its redshift is otherwise unknown. 

Analyzing long-slit spectra of both the stars and 
ionized gas, Ciri et al. (1995; see also Jore 1997) 
found that all of the gas in this galaxy is counterro- 
tating relative to the stars. Recent CO observations 
by Garcia Burillo et al. (1998) reveal the presence of 
a massive molecular nuclear disk associated with the 
ionized gas. The H I has been mapped previously by 
van Driel et al. (1989) with the Westerbork Synthesis 
Radio Telescope (WSRT). 

As evident in Figure 5, NGC 3626 is a ringed 
Sa galaxy with smooth, faint spiral structure and a 
prominent dustlane to the west of the nucleus. The 
three left panels of Figure 5 show different represen- 
tations of the R band image, displayed with different 
contrast levels and scaling to emphasize the low s\ir- 
face brightness outer disk (top), the ring (center) and 
the dustlane (bottom). Unfortunately, the telescope 
suffered poor tracking when the B and V images were 
obtained but, while the PSF is degraded, colors can 
still be derived. The B — V color image shows no 
large scale structure in the disk aside from the dust- 
lane. The reddening seen in both B — V and B — R 
at -|-10" along the minor axis coincides with the dust- 
lane. The bulge is small, dominating the inner R < 
9", with a.B — V color of ~ 0.7; two larger blue {B — V 
~ 0.4) components are found about 6" to the NW and 
SE of the bulge. As shown in the top right panel of 
Figure 5, the R band surface brightness profile is well 
fit by an n = 2 bulge, yielding a modest B/Dr ~ 
0.37. Similar values are obtained for the other bands 
as given in Table 4. Structure seen in the surface 
brightness profile coincides with the crossing of the 
dustlane and ring. 

Zooming on the central arcminute, the lower right 
panels in Figure 5 show the uniform weight H I col- 
umn density map superposed on the R band image 
(center right) and on the Ha image (bottom right). 
The grayscale in the R band image is chosen to high- 
light the dustlane west of the nucleus. As seen in 
Figure 5, Ha is detected only in the central regions, 
peaking at the nucleus and also in a ring of radius 
^15". The ring of star formation is coincident with a 
shoulder on the R-band profile between 12 and 20". 
The H I in the inner regions avoids the nucleus but 
follows closely the Ha: within the limits of the VLA 
map resolution. Comparison of Figure 5 with Fig- 
ure 1 of Garcia-Burillo et al. (1998) reveals that the 
molecular gas peaks in the same overall ring, but just 
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Table 8 
H I Properties 







N 3626 


NPMIG 


N 3900 


N 4138 


N 4772 


N 5854 




U 6341 


U 6343 


+18.0219 


U 6786 


U 7139 


U 8021 


U 9726 


Scant (mjy) 


< 0.3 


9.9±0.1 


< 0.3 


1.2+0.1 


32*' 


3.4+0.2 




< DEF >" 




+0.07 




+0.02 


+0.08 


+0.20 


1.62 


J Shi dV (Jy km s"^) 


2.0±0.2 


8.3+0.3 


1.7+0.2 


16.7+0.7 


20. 6*= 


15.1+0.6 


0.30== 


Mhi (10^ Mq) 


0.27 


1.1 


0.23 


2.9 


1.2'' 


0.9 


0.04 


Mhi/Lb (Mq/Lq) 
V21 (km s-i) 




0.06 




0.2 


0.2 


0.1 


.01 


1645 


1484 


1466 


1801 


888'' 


1040 


1654'= 


AV (km s-')'' 


139 


382 


96 


450 


340'' 


463 


128"= 


RH1/R25 




2.6 




1.8'= 


2.5'' 


1.9 




Dynamical center (1950) 


111722.8 


111726.1 


111807.3 


114633.9 


120658.3'' 


125056.0 






+183202 


+183748 


+182938 


+271757 


+435749'' 


+022621 




Inclination (°)^ 




44-68 




65 


46 25'' 


65-58 




Vmax (km s"^) 




195 




220 


208'' 


240 





°'HI deficiency measure from Magri (1994). 
''Prom Jore et al. (1996). 

°Based on reanalysis of digital archive copy of Magri's spectrum. 

''Width at 20% of peak of global profile, uncorrected for smoothing or turbulence. 

''For main disk; including the southern extension, RH1/R25 ~ 4.5 

^Where a range is given, a flat rotation curve has been assumed, and the range indicates the variation from 
the inner to the outer disk. See Jore (1997). 
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where the H I and Ha are not seen, and vice versa. 

Figure 6 shows representations of the H I distri- 
bution and velocity field extracted from the low res- 
olution (natural weight) column density map. In the 
three panels, the optical boundaries D25 x 6.25 are 
indicated by the ellipse, oriented along the optical 
major axis at PA = 157°. The left panel shows the 
H I contours and their corresponding grayscale. The 
center panel shows the H I column density contours 
superposed on a grayscale representation of the R 
band image with the same scaling as in the center 
left panel of Figure 5. The H I is found in the central 
regions at R < 60" and in a low column density ring 
of Nij ~ 1-3 xlO^o atoms cm'^ peaking at R ~ 190". 
This outer ring is both exterior to the optical radius, 
significantly more elongated, and oriented with a dif- 
ferent position angle, nearly north-south. The ratio 
of H I to optical size is large: RH1/R25 = 2.6. The 
rightmost panel shows the H I velocity field derived 
from the natural weight data cube. The velocity field 
implies fairly regular rotation in the outer ring, with 
the north side of the galaxy is approaching and the 
south, receding. Many of the contours are jagged due 
to the patchy distribution of the gas. It is also appar- 
ent that the inner H I disk has a significantly different 
kinematic position angle than that of the outer disk 
or as measured from the overall gas distribution. 

Two nearby galaxies were detected in the field of 
view of the H I map: UGC 6341, which is located 
5.'9 south of NGC 3626, and a faint object identi- 
fied as NPMIG -hl8.0219 (Klemola et al. 1987) 12.'9 
(90 kpc) to the SE. In fact, the highest H I column 
density in the field is found coincident with the faint 
companion UGC 6341, an inclined low surface bright- 
ness Magellanic spiral. The Ha image did not extend 
out this far but in the broadband images, two central 
knots are visible which makes fitting central ellipses 
impossible. There is no bulge evident. One of the 
knots appears to be the center of light (nucleus) and 
the other lies about 10" to the NW. We derive m^ 
= 15.04 ± 0.04, with e = 0.72 and PA = 153°. The 
global HI properties of these objects are summarized 
in Table 8 and the global H I profiles of all three ob- 
jects are shown in Figure 4. 

In the central regions, the uniform weight H I map 
is deep enough to show the column density contours, 
but the S/N is too low to trace the velocity field. 
Within the resolution, the dynamical center coincides 
with the nucleus, but the inner velocity field is clearly 
distorted. To the north, the velocity contours close 



at R ^12 25" but the velocity begins to rises again 
outside 45". This rise and subsequent fall appears to 
coincide with that seen in the optical emission line gas 
rotation curve, shown in Figure 7. To the south, even 
in the center, the velocity contours change continu- 
ously. Analysis of the velocity field was accomplished 
by following the common practice of fitting a tilted 
ring model, assuming circular rotation. In the center, 
the position angle of the H I is closely aligned with the 
optical axis, starting at 160° but swiftly increases to 
195° at a radius of 65". The position angle then drops 
to 175° by a radius of 140" and changes little to the 
edge of the H I disk. Figure 7 shows position velocity 
diagrams extracted from the natural weight data cube 
at position angles of 157° (upper left), along the op- 
tical major axis, and 178° (upper right), along the 
kinematic major axis of the outer H I ring. In the left 
panel, the optical stellar and emission line velocities 
derived from the long-slit spectra are superposed. It 
is clear that the Ha and H I gas show similar rotation 
curves in the inner part, but that both counterrotate 
with respect to the stellar component. Moreover, it 
is also clear that the inner rotation curves of the ion- 
ized and neutral gas are both distinct and shifted in 
orientation from that traced by the outer H I ring. 
The rotation curve of the H I in the outer ring along 
PA = 178° continues to rise to its outermost point; 
the H I velocity gradient of the outer ring is mini- 
mized along PA =90°. A large velocity gradient is 
evident within the inner ~20" in all cuts through the 
H I cube, due, at least in part, to beam smearing. 

The H I velocity contours suggest that the incli- 
nation of the H I disk varies considerably. The best 
fit inclination given by the tilted ring model starts 
at 44° and then drops to 33° out to a radius of 80". 
However, this drop is likely an artifact of the discon- 
nect between the inner H I ring and the outer H I ring; 
there is no H I in this region to constrain the veloc- 
ity field. Towards the outside, the inclination appears 
to increase continually with radius until it reaches a 
maximum of 68° at the edge of the H I disk. That 
the outer ring is more highly inclined that the stellar 
disk is also suggested by the differences in ellipticity 
seen in the isophotes of optical light and H I column 
density. With this variation in inclination, the ro- 
tational velocity remains relatively constant at 195 
km s^^ throughout the H I distribution. 

Van Driel et al. (1989) obtained an H I synthesis 
map of this object with the WSRT. The signigicantly 
lower resolution (13" x 41") and sensitivity of that 
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map did not allow discrimination of the structural 
details but there is rough correspondence to the fea- 
tures reported here. In fact, despite the poor resolu- 
tion of their derived velocity field (40" x 41"), van 
Driel et al. noted the marked shift of position angle 
between the optical major axis and the morphological 
and kinematic position angles of the outer H I ring. 

Ciri et al. (1995) obtained both absorption and 
emission line spectra at four position angles: 25, 68, 
113 and 158°. They found that the stellar com- 
ponent entirely counter-rotates with respect to the 
stars. Along the PA = 158° axis, the stellar lines 
were traced to ±30", while the ionized lines c;xtended 
to ±20". The Palomar spectra were obtained along 
PA = 157° (major axis) and PA = 67° (minor axis); 
as shown in the lower panels in Figure 7, the extracted 
rotation curves for both components are qualitatively 
similar to those shown in Figure 1 of Ciri et al. but 
extend further in the spatial dimension. The stellar 
curves rise with some deviation at about 8" to flat- 
ten at R ^20" and beyond. The stellar curves can 
be traced in these spectra to 45" on either side of 
the major axis; this radius is only slightly more than 
0.5 R25. 

Emission line rotation curves can be traced to ± 
20" in [OIII] and [SII] and to ± 30" in Ha and [Nil]. 
Although in the very center the Ha is likely affected 
by absorption, there appears to be a hint that the 
Ha emission retains its high velocity almost to the 
center, showing velocities of ~150-180 km s^^ within 
the inner 5". In fact. Ho et al. (1997a) classify 
NGC 3626 as a possible LINER, based on line ratios; 
the FWHM of the [Nil] within their window along 
PA = 170° (Ho et al. 1995) is large, 413 km s^^. The 
[Nil], [SII] and [OIII] lines show steeply rising veloc- 
ities along the major axis, but also evidence of line 
broadening or a secondary component is seen at sym- 
metric points within R < 6". In particular, the [Nil] 
lines appear to double within the innermost region 
and have been fit by two separate Gaussian compo- 
nents, as indicated in Figure 7. When fit separately, 
the second component appears to match the veloc- 
ity of the stellar component. In addition, both the 
Ha and [Nil] curves peak at symmetric radii at R ~ 
10-15" but then show a clear decline outside R ~ 
20". This decline is suggested but not evident in the 
lower sensitivity data of Ciri et al. (1995). As dis- 
cussed above, these features are also reflected in the 
H I velocity field. The small deviation from constant 
velocity seen in the stellar values along the minor axis 



is mostly likely due to a slight offset from the true 
dynamical center (i.e. pointing) or minor axis posi- 
tion angle. Along the minor axis, both the [Nil] and 
[SII] spectra exhibit double-peaking within the inner 
R ~ 6". Within the resolution limits, it is uncertain 
whether the line broadening is merely a LINER phe- 
nomenon, excess streaming, or a secondary gas com- 
ponent co-rotating with the stars. Garcia-Burillo et 
al. (1993) likewise find suggestive evidence for a sec- 
ondary gas component associated with the dominant 
stellar disk in the central regions. 

In summary, the kinematic evidence for NGC 3626 
suggests a large scale decoupling of the stars and gas. 
The H I distribution is found in two kinematically 
distinct and spatially separate rings that show a large 
difference in position angle and inclination. Both the 
morphological and kinematic position angles of the 
outer H I ring are shifted by > 20° with respect to 
starlight and the inner H I gas. Evidence of a gas 
component co-rotating with the stars is suggested, 
as hinted also in the CO map of Garcia-Burillo et 
al. (1998). 

3.2. NGC 3900 

NGC 3900 is a ringed Sa galaxy with no nearby 
companions. The nearest known object of similar red- 
shift, UGC 6791, lies to the SE at a projected sepa- 
ration of 17.1', 135 kpc at the adopted distance of 
27 Mpc. NGC 3912, also considered part of a group 
LGG 252 by Garcia (1993) lies 34.8' (275 kpc) to the 
SW. There are several small galaxies nearby; two of 
them are included in the catalog of dwarfs found by 
Binggeli et al. (1990). None of these possible compan- 
ions has a known redshift. For one, #44 of Binggeli et 
al, we failed to detect any lines in a 600 sec exposure 
obtained with the Palomar 5m plus Double Spectro- 
graph using the setup described in Section 2.2. 

As evident in the B band images shown in Figure 8, 
NGC 3900 shows a small circular bulge with a promi- 
nent blue inner ring at a radius of 34". In the exterior, 
the light from the disk abruptly decreases at roughly 
half the optical radius. Flocculent spiral structure is 
clearly evident in a diffuse, faint disk that, as seen in 
the upper panel of Figure 8, appears to be more elon- 
gated than the higher surface brightness portions, as 
might be expected for a tightly wound spiral seen at 
moderate inclination. Only a hint of Ha is evident 
in the narrow-band images with a single concentra- 
tion near the center and a V-shaped structure in the 
north corresponding to blue enhancements (B — F ~ 
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0.6) in the northern section of the ring evident also 
in the color maps. A dustlane is seen on the southern 
side looping to the east and wrapping nearly halfway 
around the galaxy. The dustlane is easily visible on 
the northern and southern ends in the B band image, 
just interior to the inner ring and more easily distin- 
guished in the cast where it overlaps the inner ring in 
the color map images. 

From examination of a HST WFPC2 F606W im- 
age, CaroUo et al. (1997) classified N 3900 among the 
early spirals with a "classical spheroidal bulge" ; the 
bulge component identified in the HST has a smaller 
effective radius than the limited resolution ground 
based images could distinguish from a central com- 
ponent. In fact, as evident in Figure 8, a n = 2 bulge 
fits quite well and gives a modest B/D/j=0.30. Struc- 
ture in the surface brightness profile is evident in all 
the broadband images where crossings of the ring and 
dustlane occur, and where the outer faint disk con- 
tributes. Use of a higher order bulge in our R band 
image yields a poorer fit with a higher B/Dfl=0.71 
in which the bulge clearly overshoots the observed 
central light profile; the n = 2 bulge also provides a 
better fit to the outer (R> 80") disk. 

NGC 3900 was mapped previously with the WSRT 

by van Driel et al. (1989). In their map, most of the 
H I is found in a ring of radius 1', with some sug- 
gestion of the western extension. This ring is clearly 
evident in Figure 9 which shows different represen- 
tations of the H I map and velocity field extracted 
from the natural weight VLA data cube. The H I in- 
deed peaks in a ring at R 50", just outside the 
prominent blue ring, and in the region dominated by 
flocculent spiral structure. The lower sensitivity of 
the WSRT map completely missed the extended dif- 
fuse structure to the north and south, evident here in 
Figure 9. In addition to the higher density H I which 
is roughly coincident with the optical galaxy, diffuse 
H I can also be traced in patches to the north and 
south. To the south, an extension is visible, clumping 
about 5' (40 kpc) from the optical center. Further- 
more, the H I column density distribution associated 
with the main body is clearly asymmetric, extending 
nearly twice as far to the west as to the east. 

The velocity contours shown in the right panel of 
Figure 9 and the top left panel of Figure 10 clearly 
suggest that the southern gas cloud is simply an ex- 
tension of the main H I distribution, the column den- 
sity dropping below the sensitivity limit of the VLA 
map in between. In the inner parts, the contours in 



the H I disk are very symmetric about the major axis 
and show regular rotation. The main H I disk extends 
only 2f9 (23 kpc) from the center; along the major 
axis, RH1/R25 = 1-8. However, comparing the total 
extent of H I gas on the southern side of the galaxy 
to the stellar disk gives RH1/R25 = 4.5; the H I in 
the south can be traced to R ~ 50 kpc. Along the 
minor axis, to the west, the H I extends about twice 
as far as to the east and deviates in velocity as would 
be expected if the H I disk were strongly warped in 
the outer parts. 

The velocity field was modeled using tilted rings 
and including the southern low column density gas 
under the assumption that it lies in the same plane, 
at a constant inclination of 65°. In the best-fit model, 
the position angle varies from 178° to 186° in the main 
H I disk and then settles back to 176° at a radius of 
200". With these parameters, the rotational velocity 
rises to -^200 km s^^ in the inner 60" and is then rel- 
atively constant except between 240" and 300" where 
the parameters are poorly constrained because there 
is little gas in that region. In these models, the west- 
ern extension is not well fit because of the asymmetry 
in the distribution along the minor axis. 

Figure 10 shows the position-velocity diagrams for 
the H I distribution along both the optical major 
(PA = 2°) and minor (PA = 92°) axes; the veloc- 
ities extracted from the long-slit spectra for stellar 
and gaseous components are superposed. The optical 
and H I rotation curves are similar, but the optical 
ones do not conclusively show the flattening seen in 
the H I data. There is also a suggestion of a spread 
in the velocity along PA =2° to the north at R ~ lf5, 
showing a blueshifted cloud with AV ~100 km s~^. 
In the upper right panel of Figure 10, the minor axis 
position-velocity diagram shows clear evidence of the 
H I ring and central hole. It also reveals a small 
change in velocity along the western H I extension, 
likely due to the proposed warp. 

As illustrated in the lower panels of Figure 10, the 
optical spectra yield both gas emission and stellar ab- 
sorption lines extending along the major axis to R 
~ ± 55". The stellar Mg I b lines are relatively weak 
outside the bulge dominated region, at radii > 10" to 
the north and > 15" to the south. The stellar compo- 
nent to the N of the nucleus shows an asymmetric line 
of sight velocity distribution that could be cither two 
separate components but is more likely simply due 
to true asymmetric deviation from a Gaussian distri- 
bution; this detail will be discussed elsewhere (Jore 
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et al. in preparation). The central Hq is strongly af- 
fected by absorption within the inner 20" to the north 
or 15" to the south, but matches nicely the [Nil] fur- 
ther out. In the inner rotation curves derived from 
the [Nil], [OIII] and [SII] lines, a significant anomaly, 
indicating blueshifts of ~160 km s~^, is seen, peak- 
ing ~5" to the north of center along the major axis. 
Beyond this radius, the velocities fall and resume the 
expected rise seen to the south at about 10". There 
is a suggestion that the [N II] line is double peaked 
in the center. Ho et al. (1997a) classify NGC 3900 
as a possible LINER based on line ratios extracted 
from a spectrum obtained along a PA of 178° (Ho et 
al. 1995). However, as can be seen in the lower panel 
of Figure 10, the [Nil] lines exhibit a systematic sinu- 
soidal wobble within ±7" from the center along the 
minor axis. While this behavior might be expected 
from a small offset in position or position angle from 
the true kinematic center, the stellar velocities, ex- 
tracted from the same spectra, remain flat. The color 
maps show a coincident feature — slightly bluer on 
the east and redder on the west of the nucleus at com- 
parable radii. 

While the main disk of NGC 3900 shows motions 
dominated by circular rotation in a quiescent pattern, 
anomalous emission line velocities imply non-circular 
motions in the central region. Interior to the optical 
edge, the H I is found in a ring or tightly wound spiral 
pattern with symmetric peaks to the north and south 
along the major axis in the region of flocculent opti- 
cal spiral structure. At low column densities, Nh < 
2 xlO^° atoms cm~^, the H I distribution is huge. 
The outer H I can be traced to a radius more than 4 
times the optical size in a diffuse and patchy distribu- 
tion that appears kinematically connected to the main 
galaxy. Along the minor axis, the H I distribution is 
both warped and strongly asymmetric. Models of the 
evolutionary history of this galaxy must account for 
the kinematic quiescence of the patchy and asymmet- 
ric H I far beyond the optical edge. 

3.3. NGC 4772 

NGC 4772 is an outlying member of the Virgo clus- 
ter, probably associated with the subclump around 
NGC 4472. The adopted distance assumes this asso- 
ciation and the derived subgroup distance of 16 Mpc 
from distance moduli given in Ferrarese et al. (2000) . 
Its nearest companion, CGCG 015-036, a compact 
high surface brightness galaxy with a systemic veloc- 
ity of V© = 844 km s'^, is found at 18'.5, 85 kpc 



at the adopted distance of 16 Mpc; to the E SE. The 
pair NGC 4809/4810 at similar velocity lies 35' (160 
kpc) to the NE. 

Dubbed the "eye galaxy," a descriptive term ade- 
quately describing its appearance as shown in Figure 
11, NGC 4772's center is dominated by a slightly flat- 
tened bulge surrounded further out by a segmented 
ring. The western and eastern sides of the ring 
form continuous structures; however, the northern 
and southern ends are broken and patchy. The inner 
galaxy shows a clear dustlane surrounding the bright 
bulge and oriented parallel to the major axis. The 
bright, red circular bulge, of B — F ~ 1.0, is clearly 
separable within the inner 5", cutoff on the north and 
east by the dustlane. Structure (filaments or braids) 
are seen within the dustlane, which must trace spiral 
structure, not simply a ring. The bulge-disk decom- 
position yields a moderate B/Dr = 0.44. 

The Ha shows a central peak, coincident with 
the nucleus, surrounded by an elongated (b/a ^0.3) 
Ha ring which itself is crisscrossed but entwined in 
the dustlane and which coincides with the H I col- 
umn density peaks. As evident in the upper panel in 
Figure 11, there are two pronounced drops in surface 
brightness, one at about 65" and a second at about 
150". This latter coincides with the outer edge of 
the inner H I ring, just about where there are sug- 
gestions of the start of spiral arms. At that radius, 
the isophotes, ^ib = 25.9, fin. = 24.0 magarcsec"^, 
become rounder, changing from e ^ 0.47 to e ~ 0.29, 
similar to what is seen in the H I distribution. At the 
faintest isophotes, this rounder low surface brightness 
distribution can be traced to more than 200", twice 
R25, and roughly to the outermost H I contours. The 
outer low surface brightness disk is blue, B — V ^ 
0.67 versus B — V ^ 0.75 - 0.85 further in, but con- 
tributes only a small amount of light (0.19 mag in B) 
to the total luminosity. The faint optical light is vis- 
ible on the DSS image but makes the galaxy rounder 
than published axial ratios suggest. It is reminiscent 
of a similar feature detected by Buta et al. (1995) in 
deep photometry of the counter-rotating Sab galaxy 
NGC 7217. 

As evident in Figure 12, the H I gas is distributed 
in two distinct rings surrounding a central hole. Co- 
incident with the Ha ring evident in Figure 11, the 
inner high column density H I ring peaks at R ^ 60", 
with a position angle of 144°, similar to that of the 
bright optical disk and with comparable eccentricity. 
Surrounding the central hole, the full ring is charac- 
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tcrizcd by moderate column densities, N/f > 7 xlO^" 
atoms cm~^. Clearly separated from this inner ring, 
the rounder outer H I ring can be traced to ~ 200" in 
radius oriented along a position angle of 150° and, as 
seen in the lower right panel of Figure 12, extends over 
an area similar to, but of different eccentricity, with 
the faint blue stellar light. The outer ring is patchy, 
but, in the north, can be traced over 90° in azimuth 
at Njj ~ 2 xlO^° atoms cm~^, with localized patches 
of 3 xlO^" atoms cm~^. The velocity field shows a 
similar discontinuity. Near the edge of the inner ring 
and in the outer ring, along the major axis, the isove- 
locity contours bend, suggesting a warp in the outer 
regions. Figure 13 shows the position-velocity dia- 
grams along the optical major (PA=147°) and minor 
(PA=57°) axes; the velocities derived from the optical 
spectra in the same regions are shown. The velocity 
rises linearly in the inner regions and then flattens out 
at a radius of 90". Within the outer ring, the H I ve- 
locities show a significant turn-over at large radii in 
the SE. Along the minor axis, the contours appear 
regular and symmetric surrounding the central hole, 
but the gradient of the outer H I velocities reinforces 
the notion that the outer gas disk is warped. 

As in the cases of NGC 3626 and NGC 3900, the 
fitting of tilted ring models is complicated by the dis- 
continuity in the H I distribution between the inner 
and outer rings. The inner region shows a rotational 
velocity smoothly rising to 240 km s~^ at a radius of 
100". If a flat rotation curve, set by the fit to the inner 
disk, is adopted, the position angle starts at 143° and 
rises steadily to 160° near a radius of 140" remain- 
ing roughly constant at -^158° further out. However, 
just at the edge of the inner H I ring, the inclination 
becomes unconstrained and discontinuous. A second 
fit was performed to smooth over the discontinuity in 
inclination by fixing the inclination for the transition 
rings to allow a smooth variation in the inclination. 
In that case, the inner rotational velocity rises much 
in the same way as in the previous model, but in the 
outer ring, it rises by 60 km s~^. Yo maintain a flat 
rotation curve, the inclination must drop outside R> 
110" from about 65° to about 58°, as suggested also 
by the rounding of the aspect ratio of the outer col- 
umn density contours. In terms of residuals, both 
of these models fit the observed velocity field equally 
well. Further details of the resultant mass modelling 
will be discussed elsewhere. 

Figure 13 also shows the stellar and ionized gas 
velocities extracted from the long slit spectra. The 



stellar lines can be traced from about 60" on the SE 
to about half as far on the NW. The ionized gas dis- 
tribution is extremely patchy. [Nil], [OIII] and [SII] 
emission is detected within the central R < ±6", and 
on the southeast side along the major axis in clumps 
at R ~ 25" and 50". No ionized emission is detected 
to the NW. The Ha line is depressed in the center 
due to absorption as far out as 15" but the outer 
Ha emission coincides nicely with that of the other 
ionized species. As also evident in the top panels of 
Figure 13, the optical rotation curves do not extend 
sufBently far to sample the rotation curve peak. 

Ho et al. (1997ab) classify NGC 4772 as a Seyfert 
Syl.9, based on line ratios and the appearance of 
broadened Ha seen in their spectrum obtained along 
PA = 5°. They also remark on the flat or double- 
peaked nature of the narrow lines. The Ha line broad- 
ening is also evident in our spectrum but the other 
species, particularly [OIII], show a systematic rise to 

80 km at R~ ± 6", more steeply and of op- 
posite sign than the stellar rise in the same region. 
It is important to note that, unlike Ho et at, who 
chose to orient the slit close to the parallactic an- 
gle, our setup placed the slit along the optical axis 
at P.A. = 147°. Furthermore, as evident in Figure 
13, the minor axis [Nil], [OIII] and [SII] velocities 
rise rapidly in the inner 5" up to a velocity of 100 
km s^^ and then decrease out to 15" in radius. In 
this case, the [OIII] provides the best tracer, showing 
a velocity peak along the minor axis at 4- 5" followed 
by a decline back to the systemic by the time it dis- 
appears at about 15". At the same time, the minor 
axis stellar velocities are very flat (verifying that the 
slit was correctly aligned along the minor axis of the 
stellar component). The kinematic signature suggests 
the presence of a misaligned disk or bar. The only 
unique feature seen in the color maps is a flattening 
of the B — R in the central region along the minor axis 
and a hint of smaller plateaus in all the color maps 
along the major axis. These two color features occur 
on comparable scales as the central rotation curve pe- 
culiarities. 

NGC 4772 may represent the end state of a pro- 
grade merger in which the transfer of angular momen- 
tum leads to an outward spread of the disk (Quinn et 
al. 1993). The steep rise of the ionized gas velocities 
along the minor axis and the apparent decoupling on 
gas emission along the major axis suggests the pres- 
ence of a misaligned embedded disk or bar. Reminis- 
cent of NGC 3626, the H I is contained in two separate 
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but concentric rings of differing position angle and ec- 
centicity. The inner ring coincides witli the sites of 
current star formation in the optically bright inner 
disk, while the other rounder one, wholly outside R25 
but coincident with a faint blue stellar component. 

3.4. NGC 5854 

NGC 5854 is a gas poor Sa with no evidence of cur- 
rent star formation. Relatively isolated, it is likely an 
outlying member of the NGC 5846 group (Haynes & 
GiovaneUi 1991; Garcia 1993; Zabludoff & Mulchaey 
1998). Its nearest neighbor is CGCG 021-009 at 33.'1, 
230 kpc at the assumed distance of 24 Mpc. An il- 
lustrative summary of results is presented in Figure 
14. 

The B band image shows an inclined galaxy with a 
large elongated bulge and little apparent spiral struc- 
ture. Some hint of arms is seen in the two small 
stumps on the NE and SW edges of the bulge. Un- 
fortunately, the B band image was taken under non- 
photometric conditions so that the color information 
is limited to V — R. In the V — R image, little struc- 
ture is seen, although the galaxy is clearly redder 
than the others. It was previously imaged by Bal- 
cells & Peletier (1994) who find B/Br = 0.46 using 
a 1-D decomposition method and find similar flatness 
of the color profile. As discussed below, no Ha is 
seen in the long slit spectra, but Ha absorption is 
seen across the disk and can be traced almost as far 
as the Mg I b lines. The galaxy is thus not currently 
forming stars, but its significant population of A stars 
suggests a post-starburst state. Of the galaxies dis- 
cussed here, it has the highest B/D/j ~ 0.6 for n = 2, 
but it is clear that the n = 4 fit, illustrated in Fig- 
ure 14, gives a more appropriate B/D/j ^ 0.8. No 
Ha image was obtained. 

As reported by Magri (1994), the H I mass of 
NGC 5854 is very low (- 3 x 10^ M©) and hence 
it was not feasible to map its H I with the VLA. The 
optical spectra are dominated by stellar absorption 
lines which can be traced, by both Mg I b and Ha, 
out to ±45", further than, but qualitatively similar 
to the rotation curve displayed in Simien & Prugniel 
(1997). Those authors note V„ia.i = 128±8kms-i for 
the stellar component in their spectrum traced to ~ 
27'.'5. We find a comparable value at that radius, but 
a larger Yrnax = 167 ± 10 km s"^ at R -40". The 
Mg I b rotation velocities rise more shallowly than the 
Ha lines do, possibly because of extinction. The [Nil] 
and [OIII] emission reveals a distinct counter-rotating 



gas component traceable from the center out to 6" on 
the NE side and 8" on the SW side. The maximum 
rotational velocity of this [Nil] feature is 64 km s~^ at 
~4" NE of the optical center of Hght, the [OIII] ap- 
pears in a ring coextensive with the [Nil] gas but the 
velocities are higher, reaching ~ 100 km s~^. The 
stellar velocity profiles along the minor axis show a 
slight trend that most likely arises because of a small 
pointing offset. The weak [Nil] shows some disper- 
sion but is essentially flat, within the errors; no other 
ionized emission is seen. 

Magri (1994) detected a weak and narrow H I sig- 
nal from NGC 5854 at Arecibo after more than 2 
hours of on-source integration. A zoomed reproduc- 
tion of Magri's spectrum extracted from our digital 
archive is shown in Figure 14. The systemic veloc- 
ity of the H I is 1654 km s^^ which matches the 
value of 1663 ±10 km s~^ given by Simien fc Prug- 
niel (1997). We have remeasured the H I line width 
in the digital spectrum and corrected it for smooth- 
ing, etc. following Haynes et al. (1999b) to obtain a 
value, uncorrected for inclination and turbulence, of 
W21 ~95 km s~^. This value is too small compared to 
the full rotation width of the galaxy, but is compara- 
ble to the velocity spread seen in the [Nil] in the inner 
disk. A simple interpretation of the narrowness of the 
H I profile is that that the H I and [Nil] both arise in 
a small central disk of R < 8" that counter-rotates 
relative to the stars. At the same time, the primary 
(stellar) disk is seen in a post-starburst phase, show- 
ing strong Balmer absorption throughout. The cur- 
rent gas-poor state of NGC 5854 may thus be result 
of the accretion of a relatively dense satellite that has 
heated the disk and smoothed the spiral structure, 
leaving only the counter-rotating gas core as evidence 
of its occurrence. 

3.5. Another Look at NGC 4138 

As mentioned previously, the first object mapped 
in H I with the VLA as part of this study was 
NGC 4138. A discussion of the optical and H I spec- 
tral observations is presented in Jore et al. (1996); the 
gas counterrotates with respect to the primary stel- 
lar component, but a secondary stellar popiflation, 
kinematically coupled to the gas, is visible. The am- 
plitude of the crosscorrelation function of the stellar 
lines coincides spatially with an Ha ring, previously 
noticed by Pogge & Eskridge (1987), at a radius of 

20". The two stellar components are characterized 
by distinct velocity dispersions, implying separate ori- 
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gins. As part of the current study, we added _B, V, R 
and Ha images to the dataset for this morphologically 
normal but kinematically unusual object. 

As evident in the top left panel of Figure 15, 
NGC 4138 is a small bulged, smooth armed Sa galaxy 
with a prominent dustlane in the SE quadrant. A 
foreground star is located at a radial distance of 
4" from the center along the major axis. As seen 
in the middle panels of Figure 15, the inner H I forms 
a partial ring, coincident with the Ha ring, on the 
north and east, surrounding a central hole. In. B — V, 
the galaxy shows a distinctly red bulge {B — V ^ 
1.1) and a blue {B — V ~ 0.6) ring located in the 
approximate position of the Ha ring. The blue ring 
is particularly strong on the northern and southern 
portions of the galaxy with breaks in the ring on the 
east and west. The eastern break is due, at least in 
part, to the dustlane. The V — R image shows sim- 
ilar structure. As shown in the upper right panel of 
Figure 15, the bulge-disk decomposition gives a small 
value: BfDR = 0.17. 

There are several outstanding features in the color 
profiles which correlate well with features in its rota- 
tion curve as presented in Jore et al. (1996). At ±13", 
where the secondary component of stars begins, until 
±35" when the component is no longer detectable in 
the long slit spectrum, there is an obvious reddening 
in the major axis profiles. There is an abrupt tran- 
sition moving from the center into the ring on the 
southeast side and an abrupt transition out of the 
ring on the northwest side, There are no apparent 
variations in the color profiles along the minor axis 
aside from a peak redder color in the center. 

The relationship of the H I gas reported previ- 
ously by Jore et al. to the ring and the dustlane is 
demonstrated here in more detail in Figure 15. In- 
deed, the relationship of the dustlane and Ha ring to 
the H I distribution seen in NGC 4138 is very similar 
to that seen in NGC 3626 (Figure 15 versus Figure 5), 
the other case of large-scale counter-rotation. Like- 
wise, the H I extent is very large, Rhi /R25 = 2.5 and, 
as discussed in Jore et al., the velocity contours along 
the minor axis suggest a strong warp beyond R25. In 
contrast to the H I distributions seen in NGC 3626 
and NGC 4772, however, the outer HI does not form 
a distinct and coherent ring, but rather is extremely 
patchy beyond the optical edge. 

NGC 4138 remains the sole example in our Sa sam- 
ple of a galaxy with an extended counter-rotating 
stellar component. Nonetheless, it bears similarities 



to the other objects in terms of its morphological 
"boring-ness" , relative isolation, large, low surface 
density H I extent, and overall regularity of its ve- 
locity field. 

4. Discussion 

In this paper, we have presented a combined imag- 
ing and spectroscopic dataset for four Sa galaxies; im- 
ages are also included for a fifth, NGC 4138, discussed 
previously by .Jore et al. (1996). All of the galax- 
ies show evidence of departure from kinematic nor- 
malcy ranging from large-scale counter-rotation to 
decoupled central gas and/or stellar components. De- 
spite the difference of detail, several unifying themes 
seem critical: (1) All of the objects are relatively iso- 
lated, morphologically normal, unbarred Sa galaxies. 
(2) Rings are prominent features both in the optical 
broadband light and in the distributions of ionized 
and neutral gas. The apparent rings may also result 
from tightly wound spiral arms viewed at moderate 
inclination. (3) The kinematically decoupled gas com- 
ponents appear to be associated with sites of current 
or at least recent star formation. In all cases where 
Ha emission is detectable in the region of interest, 
it coincides well with the location of kinematically 
decoupled components. When there is no Ha emis- 
sion in these regions, the decoupled components often 
show some form of a color change, often towards the 
red. (4) The gas-rich galaxies (all except NGC 5854) 
contain moderate H I masses, but, because the gas is 
spread over a large area Rff//R25 > 2, the globally 
averaged H I surface densities < anj >= Mhi/ttR^j 
are very low, typically of order 0.5-1.0 Mq pc~^. In 
NGC 3626 and NGC 4772, the H I is found in two con- 
centric but distinct rings. In NGC 3900, the gas ex- 
terior to the optical edge is patchy, but can be traced 
50 kpc to the south, giving ~ 4.5. Along the 

minor axis, the exterior H I is strongly asymmetric. 
(5) In all four H I maps, including that of NGC 4138, 
the velocity field is dominated by circular rotation, 
but significant departures from motion in a quiescent 
disk are evident. In all, the outer H I velocity con- 
tours suggest significant warping of the outer disk, 
roughly beginning at the optical edge and following 
Briggs' (1990) rule for warps, as commonly seen in 
other galaxies. The similiarities found in the five Sa 
galaxies are consistent with the scenario that the kine- 
matically distinct gas components arise from slow, 
minor mergers of gas -rich satellites with an already- 
formed disk galaxy. Differences in their appearance 
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thus arise from differences in the circumstances of the 
merger events. Retrograde primordial gas infall, as 
discussed by Thakar et al. (1997) is also offered as a 
possible solution to produce counter rotating compo- 
nents, although seemingly more difficult to justify in 
the growing number of known cases. 

Though most numerical work on the merger pro- 
cess has focussed on equal mass progenitors, recent 
studies have begun to address the minor merger phe- 
nomenon (Quinn et al. 1993; Hernquist & Mihos 1995; 
Walker et al. 1996; Thakar et al. 1997; Bekki 1998; 
Thakar & Ryden 1998; Taniguchi 1999). Although 
much of this work is preliminary, requiring further re- 
finement of the range of galaxy characteristics and or- 
bital parameters and the influence of gas dissipation, 
triaxiality and dark matter c;ontent. wc can nonethe- 
less examine the main characteristics of the five galax- 
ies under the assumption that we may be seeing the 
results of varying acquisition scenarios. 

The extensive H I distributions seen in four of the 
five galaxies discussed here may result from the slow 
accretion of a gas-rich satellite that is tidally stripped 
before the merger occurs. Accretion of a low den- 
sity gas-rich satellite would minimize the disk heating 
problem and might favor the formation of extended 
H I disks or rings, but perhaps not inner rings. The 
formation of rings is seen in a variety of infalling gas 
models, both due to orbit crowding and to true self- 
gravitating rings. Also, Quinn et al. (1993) suggest 
that transport of angular momentum outwards will 
tend to cause the disk to spread in radius. The newly- 
acquired gas may reach sufficient densities, if clumpy, 
to form a new generation of stars. This possibility 
is suggested by the faint, blue, low-surface brightness 
outer disk seen in NGC 4772 and a similar excess of 
light found around the counter-rotating Sab galaxy 
NGC 7217 (Buta et al. 1995). 

Clearly the timcscales of the merger events respon- 
sible for the individual characteristics noted here are 
likely to vary greatly. Relatively gas-poor and large- 
bulged, NGC 5854 may represent the result of the 
accretion of a more massive satellite travelling in a 
prograde orbit (Walker et al. 1996). Merger events 
might lead to bulge build- up and even induce large- 
scale star formation, leaving the galaxy in the post- 
starburst state in which NGC 5854 is found. Further- 
more, one of the primary arguments made against 
the tidal forcing of warps is that the modes would 
damp out after a few galactic rotations. Perhaps in 
the cases that suggest strong warping of the H I, the 



minor mergers may be relatively recent events. Over 
time, the warps may damp, or be reexcited by ad- 
ditional interactions. If there are no impediments to 
star formation in a coimtcr- rotating gas component, 
NGC 4138 with its two counter-rotating extended 
stellar disks may represented the future evolutionary 
state of NGC 3626. 

The similarities among the results for these five 
galaxies lead to several questions and conundrums. 
(1) Is it key that the sample is composed of Sa galax- 
ies? The accretion of moderate amounts of gas may 
also sweep up the gas contained within the primary 
disk itself, ultimately leading to a smoothing of the 
spiral structure. Do minor mergers lead to bulge 
build-up and disk heating, driving the post-merger 
morphology towards earlier spiral types? (2) Damp- 
ing the process of disk heating is a challenge for mod- 
els of mergers, even minor ones, if the accreted satel- 
lite is sufficiently massive or dense. The resultant 
disk heating will be relatively lessened if the satellite 
is predominantly gaseous. What are the conditions 
necessary to avoid over-heating the disk and disturb- 
ing the overall velocity field and spiral structure? (3) 
Is the fact that these galaxies are all unbarred im- 
portant? Loss of angular momentum may drive the 
material inward toward the central regions where it 
accumulates in a nuclear disk as is often discussed for 
the kinematically decoupled cores seen in ellipticals. 
In the early stages of galaxy evolution, the formation 
of bars seems inevitable, but if the central concen- 
tration of mass is sufficient, the development of an 
inner Lindblad resonance will inhibit the central flow 
of gas (Sellwood & Moore 1999). Thus, since the pro- 
genitor primary was already a well-developed object, 
infalling material may be stopped at the inner Lind- 
blad resonance, allowing the build-up of gas in a ring 
with subsequent star formation. Does this process 
produce the inner Ha and H I rings seen here? Jung- 
wiert & Palous (1996) propose that unbarred ringed 
galaxies must possess weak central oval distortions, 
while the two-stream instability picture of Lovelace 
et al. (1997) for galaxies with counter-rotating com- 
ponents suggests the generation of strong m=l spiral 
waves. Are the inner kinematic peculiarities evidence 
of such oval distortions or m=l streaming? (4) Al- 
though these Sa galaxies are relatively isolated, they 
are all members of loose groups. In fact, Sa galaxies 
are found typically in denser environments than their 
later type counterparts. Is the relative morphological 
segregation seen across the spiral sequence related to 
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the likelihood of minor mergers and their increased 
probability in environments characterized by moder- 
ate density and low velocity dispersion? 

While examination of the rate of strongly disturbed 
galaxies may provide a reliable estimate of the major 
merger rate (e.g. Keel & Wu, 1995), minor merg- 
ers may be morphologically hidden, revealed only by 
the presence of their resultant kinematic peculiari- 
ties. Because of the ambiguities involved in disentan- 
gling the kinematic clues, it is not possible to trace 
with any degree of precision the evolutionary history 
of the multiple stellar and gaseous components evi- 
dent in these galaxies. Food for thought, however, is 
the possibility that, over the course of its lifetime, a 
disk galaxy might undergo multiple minor accretion 
events with kinematic as well as morphological con- 
sequences. In particular, more complicated histories 
involving multiple minor mergers are likely to produce 
heterogenous properties, particularly those related to 
properties of "form" . We suggest that the heterogene- 
ity of the Sa class in particular is the outcome of such 
complicated life histories, the memory of which is sig- 
naled only by kinematic clues such as those identified 
here. 
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Fig. 1. — Cleaned continuum-subtracted H I line 
channel maps of NGC 3626 derived from the natu- 
ral weight (low resolution) data cube. Contour levels 
range from 0.46 to 4.6 mJy beam~^ in increments of 
0.46 mJy beam'^ The beam size is 16'.'8 x 16'.'2. 



Fig. 2. — Cleaned continuum-subtracted H I line 
channel maps of NGC 3900 derived from the natu- 
ral weight (low resolution) data cube. Contour levels 
range from 0.42 to 6.42 mJy beam^^ in increments of 
1.0 mJy beam-^ The beam size is 22'.'5 x 20'.'0. 



Fig. 3. — Cleaned continuum-subtracted H I line 
channel maps of NGC 4772 derived from the natu- 
ral weight (low resolution) data cube. Contour levels 
range from 0.62 to 8.62 mJy beam~^ in increments of 
1.0 mJy beam-^ The beam size is 31'.'7 x 23'.'1. 



Fig. 6. — H I mapping results obtained from the low 
resolution cube for NGC 3626 and UGC 6341. Left 
panel: H I column density contours and grayscale 
ranging from 1.0 {2a) to 17.5x10^" atoms cm^^ in 
steps of 1.5x10^° atoms cm~^. Middle panel: The 
same H I column density map overlayed on the R 
band imago. The optical greyscale has the same shad- 
ing as in the center left panel of Figure 5. Right panel: 
The H I velocity field map. Velocity contours range 
from 1275 (light) to 1675 (dark) km s^^ in increments 
of 25 km s~^. In all three panels, the ellipse outlines 
the optical galaxy as given by £'25 x ^25 oriented along 
the optical major axis, PA = 157°. 



Fig. 4. — Global profiles derived from VLA mapping 
for the target galaxies and neighboring companions. 
H I parameters are given in Table 8. 



Fig. 5. — Imaging results for NGC 3626. Top left 
panel: R band image with the scaling set to showing 
the maximum extent of starlight. Middle left panel: 
Expanded R band view illustrating the faint ring and 
spiral structure. The ellipse has the dimensions of D25 
X d25 and is oriented along the optical major axis. 
Lower left panel: The inner 1' shown to emphasize 
the dustlane west of the nucleus. Upper right panel: 
Results of the two-dimensional bulge-disk decompo- 
sition for a n = 2 bulge. The solid line traces the 
light distribution while the disk and bulge are shown 
as long dashed and short dashed lines respectively. 
Middle right panel: The high resolution H I column 
density map superposed on a greyscale representation 
of the central R band image The first H I contour 
is at 1.0x10^° atoms cm~^. The next six contours 
range from 4.0 to 11.5x10^*^ atoms cm~^ in steps of 
1.5x10^0 atoms cm and the remainder range from 
13.5 to 18.5x10^0 atoms cm'^ in steps of 1.0x10^° 
atoms cm~^. Lower right panel: Identical contours 
of the high resolution H I column density map super- 
posed on the Ha image. 



Fig. 7. — Position-velocity cuts obtained from the 
low resolution H I cube and the optical spectra for 
NGC 3626. Upper panels: H I position-velocity di- 
agrams along four slices through the low resolution 
data cube. The top left panel shows the slice along the 
optical major axis. PA = 157°, along which the opti- 
cal spectra were obtained. In addition to the H I con- 
tours, the velocities of the optical stellar absorption 
line (open squares) and emission line (filled circles) 
are also shown. The secondary gas component is rep- 
resented by open circles. The top right panel shows 
the slice at PA = 178°, the major axis of the outer 
H I distribution. The center loft panel shows the slice 
along the optical minor axis, PA — 67°; the velocities 
of the stellar and ionized gas components obtained 
from the long slit spectra are also shown. The cen- 
ter right panel shows the slice along a position angle 
of 90°, the minor axis of the outer H I ring where 
the velocity gradient is minimized. The lower panels 
show the optical rotation curves derived for NGC 3626 
along the major axis (PA = 157°) and along the mi- 
nor axis (PA = 67°). Separate symbols show the ro- 
tation curves derived from the stellar absorption and 
gas components. Error bars are shown. In the inner- 
most regions, the [Nil] lines appear double. 
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Fig. 8— Results of optical imaging for NGC 3900. 
Upper panel: B band image with scale set to empha- 
size outer structure. The ellipse has the dimensions 
of D25 X d25 and is oriented along the optical major 
axis, PA = 2°. Middle panel: Low contrast B band 
grayscale showing spiral structure and dustlanes in 
the inner galaxy. The spatial scale is zoomed by a fac- 
tor of two relative to the upper panel. Lower panel: 
Results of the two-dimensional bulge-disk decompo- 
sition for a n = 2 bulge. The solid line traces the 
light distribution while the disk and bulge are shown 
as long-dashed and short-dashed lines respectively. 



Fig. 9. — H I mapping results obtained from the low 
resolution cube for NGC 3900. Left panel: H I column 
density contours ranging from 0.3 {2a) to 12.3x10^° 
atoms cm~^ in steps of 2.0x10^° atoms cm~^. Mid- 
dle panel: Similar H I column density map, with a 
stepsize of 1.0x10^° atoms cm~^ overlayed on the B 
band image. Right panel: The H I velocity field map 
showing velocity contours ranging from 1575 (light) 
to 2050 (dark) km s~^ in increments of 25 km s~^. In 
all three panels, the ellipse outlines the optical galaxy 
as given by D25 x ^25 oriented along the optical major 
axis PA = 2°. 



Fig. 11.- Imaging results for NGC 4772. Top left 
panel: B band image scaled to show the faint outer 
disk. Middle left panel: B-band image of the same 
region but with the grayscale set to emphasize the 
inner regions, showing that the dustlane cuts across 
the bulge in the NE. Lower left panel: Ha image. The 
ellipse superposed on the three images has the dimen- 
sions of D25 X d25 and is oriented along the optical 
major axis, PA = 147°. Upper right panel: Results 
of the two-dimensional bulge-disk decomposition for 
a n = 2 bulge. The solid line traces the light distri- 
bution while the disk and bulge are shown as long- 
dashed and short-dashed lines respectively. Middle 
right panel: The high resolution H I column density 
map superposed on a greyscale representation of the 
B band image in the inner regions. Contours extend 
from 0.6x10^° to 14x10^° atoms cm~^ in steps of 
4.0 xlO^" atoms cm . Lower right panel: Identical 
contours of the high resolution H I column density 
map superposed on the Ha image. 



Fig. 10. — Position-velocity cuts obtained from the 
low resolution H I cube and the optical spectra for 
NGC 3900. Top panels: H I position velocity di- 
agrams along two slices through the low resolution 
data cube. The top left panel shows the slice along 
the optical major axis, PA = 2°, along which the opti- 
cal spectra were obtained. In addition to the H I con- 
tours, the rotation curves traced by the optical stellar 
absorption lines (open squares) and the gas emission 
lines (filled circles) are also shown. The top right 
panel shows the slice at PA = 92°, along the mi- 
nor axis. The lower panels show the optical rotation 
curves derived for NGC 3900 along the major (PA = 
2°) and minor (PA = 92°) axes. Separate symbols 
show the rotation curves derived from the stellar ab- 
sorption and gas components. Error bars are shown. 



Fig. 12. Results for the large-scale H I mapping of 
NGC 4772. Upper left panel: H I column density con- 
tours superposed on the B band image. Contours are 
shown at levels of 0.5, 1.5, 2.5, 4.0, 6.0, 8.0, 10.0 and 
12.0 and 14.0 xl0^° atoms cm~^. Upper left panel: 
Greyscale and contours of the H I column density as 
shown in the left panel. Lower left panel: The out- 
ermost HI column density contour at 5x10^^ atoms 
cm~^ superposed on a grayscale representation of the 
B band image scaled to enhance the faint outer disk. 
Lower right panel: The velocity field of NGC 4772, 
shown as contours and grayscale extending from 770 
km s-i (NW) to 1250 km s'^ (SE) in steps of 30 
km s~^. In all panels, the ellipse has dimensions D25 
X d25 and is oriented along the optical major axis, 
PA = 147°. 
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Fig. 13. — Position velocity cuts obtained from the 
low resolution H I cube and the optical spectra for 
NGC 4772. Top panels: H I position-velocity di- 
agrams along two slices throiigh the low resolution 
data cube. The top right panel shows the slice along 
the optical major axis, PA = 147°, along which the 
optical spectra were obtained. In addition to the 
H I contours, the velocities derived from the optical 
stellar absorption lines (open squares) and ionized gas 
emission lines (filled circles) are also shown. The top 
right panel shows the slice at PA = 57°, along the 
minor axis. The lower panels show the optical rota- 
tion curves derived for NGC 4772 along the major 
axis (PA = 147°) and along the minor axis (PA = 
57°). Separate symbols show the rotation curves de- 
rived from the stellar absorption and gas components. 
Error bars are shown. 



Fig. 14.— Results for NGC 5854. Upper left 
panel: R band image. The ellipse has dimensions 
D25 X das and is oriented along the optical ma- 
jor axis, PA = 55°. Upper right panel: Results of 
the two-dimensional bulge-disk decomposition for a 
n = 4 bulge. The solid line traces the light distri- 
bution while the disk and bulge are shown as long- 
dashed and short dashed lines respectively. Center 
right panel: Central channels of the H I single dish 
spectrum obtained at Arecibo by Magri (1994). The 
lower panels show the optic;al rotation cnu'ves derived 
for NGC 5854 along the major axis (PA = 55°) and 
along the minor axis (PA = 145°). Separate symbols 
show the rotation curves derived from the stellar ab- 
sorption and gas components. Error bars are shown. 



Fig. 15.— New imaging results for NGC 4138. Top 

left panel: B band image scaled to illustrate the in- 
ner bulge, disk and dustlane. The ellipse has dimen- 
sions D25 X d25 and is oriented along the optical ma- 
jor axis, PA = 150°. Top right panel: Results of 
the two-dimensional bulge-disk decomposition for a 
n = 2 bulge. The solid line traces the light distri- 
bution while the disk and bulge are shown as long - 
dashed and short-dashed lines respectively. Center 
left panel: The inner H I contours superposed on a 
zoomed B band image. The contours arc at levels of 
20 to 95 X 10^^ in steps of 15 x 10^^ atoms cm~^. Cen- 
ter right panel: Similar display of the H I superposed 
on the H alpha image. Bottom left panel: H I column 
density contours superposed on the V band image. 
The H I map shows contours at levels of 3.4, 6.9 and 
13.8 xlO^'' atoms cm~^ plus contours from 20 to 95 
xlO^^ atoms cm~^ in steps of 15 atoms cm~^. For 
further details, see Jore et al. (1996). The ellipse is 
the same as in the upper left panel. Bottom right: 
Identical H I column density contours as the bottom 
left panel, here also shown as grayscale to emphasize 
the patchiness of the outer H I distribution. 
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